Gravitational Waves from Chaotic Dynamical System 
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To investigate how chaos affects gravitational waves, we study the gravitational waves from a 
spinning test particle moving around a Kerr black hole, which is a typical chaotic system. To 
compare the result with those in non-chaotic dynamical system, we also analyze a spinless test 
particle, which orbit can be complicated in the Kerr back ground although the system is integrable. 
We estimate the emitted gravitational waves by the multipole expansion of a gravitational field. We 
find a striking difference in the energy spectra of the gravitational waves. The spectrum for a chaotic 
orbit of a spinning particle, contains various frequencies, while some characteristic frequencies appear 
in the case of a spinless particle. 

PACS numbers: 04.30.-w,95.10.Fh,04.20.-q 



I. INTRODUCTION 



The current ground-base gravitational wave detectors 
such as LIGO, TAMA300, GEO600 and VIRGO 0,1111 
0] are now starting the "science run" and the laser space 
antenna for gravitational waves (LISA) (5| may operate in 
the near future. Gravitational waves will bring us various 
new information about relativistic astrophysical objects. 
If we will detect gravitational waves and compare them 
with theoretical templates, we may be able to determine 
a variety of astrophysical parameters of the sources such 
as their direction, distance, masses, spin, and so on. The 
direct observation of gravitational waves could resolve 
strong-gravitational phenomena such as a black hole for- 
mation. Furthermore, we may be able not only to verify 
the theory of gravity but also to find new information 
or to recover new physics at high density or high energy 
region. Therefore we need to make theoretical templates 
of gravitational waves from various astrophysical objects 
and phenomena to extract useful information from grav- 
itational waves. 

There are a lot of astrophysical objects and phenomena 
as gravitational wave sources, e.g., a closed binary sys- 
tem or a supernova explosion. Many attempts to study 
their templates have been done so far. In this paper, 
we particularly focus on a gravitational waves form a 
chaotic dynamical system. It is motivated as follows. 
Chaos appears universally in nature and it is expected 
to explain various nonlinear phenomena. Because chaos 
indicates unpredictability of motion of objects, it comes 
into a question whether or not that it is possible to ex- 
tract useful information from a chaotic system. Many 
studies have been made on chaos in mechanics and as- 
trophysics since the research on the three-body problem 



* e-mail address : kiuchi@gravity. phys. waseda. ac.jp 
t e-mail address : maeda@gravity. phys .waseda. ac.jp 



by Poncare. It is also important to study chaos in general 
relativity because the Einstein equations are nonlinear. 
Many authors have reported chaos in general relativity 
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|23 |. As for a realistic astrophysical object, there has 
been a discussion whether or not chaos occurs in a com- 
pact binary system [2ll l23l | . It seems that there is 
not enough time in a compact binary for chaos to occur 
even in a non-integrable system. Then we may not have 
so far any concrete example of general relativistic chaos 
for a realistic situation. However, we believe that it is 
necessary to study chaos in general relativity. Because 
chaos will occur in a realistic strong gravitational field 
described by general relativity. 

In this paper, we study some observational feature of 
chaos in general relativistic dynamical system. We focus 
on gravitational waves from a chaotic dynamical system. 
As a model of a chaotic dynamical system, we adopt a 
spinning test particle around a rotating black hole. It 
mimics a model of compact objects orbiting around a su- 
per massive black hole. Recent observations support the 
existence of such massive black holes and LISA will be 
sensitive to the gravitational wave from such systems. As 
for a background spacetime, we assume the Kerr metric 
since such a supermassive black hole is usually rotating. 
The previous works 0,1111 revealed that orbits of a spin- 
ning test particle around a black hole can be chaotic due 
to the spin-orbit coupling. We evaluate the gravitational 
waves emitted from this system by using multipole ex- 
pansion of a gravitational field. |2J, |2l |2l| . 

Carter showed that the equations of motion for 
a spinless particle are integrable. After his work, John- 
ston discovered that motion for some initial parameters, 
the particle motion can be very complicated as if chaos 
occurred, although the system is integrable [2l|. Then 
we also analyze such non-chaotic but complicated orbital 
motion to clarify whether or not chaos is essential in 
our results. We compare the wave forms and the energy 
spectra of the gravitational waves from both chaotic and 
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nonchaotic systems and analyze the effects of chaos on 
gravitational waves. 

This paper is organized as follows. In Sec. II we shall 
briefly review the basic equations of both a spinless par- 
ticle and a spinning particle in relativistic spacetime, a 
supplementary condition, and some constants of motion. 
Specifying the background spacetime to be the Kerr met- 
ric, we show our results in Sec. III. Summary and some 
remarks follow in Sec. IV. 

Throughout this paper we use units c = G — 1. Our no- 
tation including the signature of the metric follows that 
of Misner Thorne and Wheeler(MTW) 



II. BASIC EQUATIONS FOR A TEST 
PARTICLE 



First we explain our system more precisely. 



A. Background spacetime 

We consider the Kerr metric as a background space- 
time. In the Boyer-Lindquest coordinates, it is given by 



, , / 2Mr\ , 2 4Afarsin^6' , , , 
ds^ = - I 1 dt^ dtd(j) 



A 



-dr^ + Sdr sin^ 6* r 



2 A/aV sin^ 6 



where 



E = + cos^ 6*, 
A = - 2Afr + 



(2.1) 



(2.2) 
(2.3) 



and M and a are the mass of a black hole and its angular 
momentum, respectively. 



B. A spinless test particle 

In this subsection, we summarize the equations of mo- 
tion for a spinless test particle. It is well known that a 
spinless test particle moves along geodesic 



dT^ 



„ dx" dxP 

= 

dT dT ' 



(2.4) 



where r is an afhne parameter of the orbit. In a station- 
ary axisymmetric spacetime, a spinless particle has two 
constants of motion; the Energy E and the ^-component 
of the angular momentum L^- The particle's rest mass fi 
is also constant. Carter discovered third constant of mo- 
tion and showed that the system is integrable [23 ■ Eq. 



(|2.4() is reduced to a set of the differential equations as 

(2.5) 
(2.6) 
(2.7) 

-P, (2.8) 



d0 








^d^ 


dr 












= - {aE 






dt 


= —a{aE 





A 



where 
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C~cos^6\a^{l ~ E^) 
E{r'^ +a 



p2 _ A{r^ + (L, - aEf + C}, 



(2.9) 

(2.10) 
(2.11) 



and C is a constant discovered by Carter and called the 
Carter constant. Note that because of the presence of 
the Carter constant, the orbits of a particle will never be 
chaotic. 



C. A spinning test particle 

The equations of motion for a spinning test particle 
in a relativistic spacetime were first derived by Papa- 
petrou ,30.j and then reformulated by Dixon [Sj. They 
used the pole-dipole approximation, in which multipole 
moments of an object higher than a mass monopole and 
a spin dipole are ignored in the limit of a point particle. 
The set of equations is given as 



dr 2 



dr 



2p^Pv 



(2.12) 
(2.13) 
(2.14) 



where t, , and S^'^ are an affine parameter of the 
orbit, the four-velocity of a particle, the momentum, and 
the spin tensor, respectively, deviates from a geodesic 
due to the coupling of the Riemann tensor with the spin 
tensor. We also need an additional condition which gives 
a relation between and v^. We adopt the condition 
formulated by Dixon [sj 



P^.S^" = 0. 



(2.15) 



This condition consistently determines the cetnter of 
mass in the present system. Using Eq ()2.15f) . we can 
determine the relation between p^" and v^" explicitly as 
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where 



(2.16) 
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and iV is a normalization constant, which is fixed by a 
choice of the affine parameter r as we will explicitly show 
it later, = /n is a unit vector parallel to the four- 
momentum p^, where the mass of the particle jj, is defined 

by 



(2.18) 



To make clear the freedom of this system, we have to 
check the conserved quantities. Regardless of the sym- 
metry of the background spacetime, it is easy to show 
that fi and the magnitude of spin S defined by 



q2 ^ Q CM^ 



(2.19) 



We fix the afhne parameter r using the condition w^m^ = 
~S. This gives N as 

N = l + ^*R*SyPS'^p\ (2.29) 



Our conserved quantities are also rewritten as 

-2 - ^ CM (2.30) 
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Q = - ^^^■,.e'^'"'^PpS.. (2.31) 



III. NUMERICAL RESULTS 



are constants of motion. If a background spacetime pos- 
sesses some symmetry described by a Killing vector 



(2.20) 



is also conserved |31|. For a spacetime with both axial 
and timelike Killing vectors such as Kerr spacetime, we 
have two conserved quantities, i.e., the Energy E and 
the z-component of the total angular momentum of a 
spinning particle. 

For numerical reasons, we rewrite our basic equations 
as follows. Following the appendix in jl^ . we write down 
the equations by use of a spin vector 



(2.21) 



where e^"''"' is the Levi-Civita tensor. The basic equa- 
tions are now 



dxP 

Dpi' 1 



(2.22) 
(2.23) 



ClT fl 

_^ = pM^*^ ^-^P^V, (2.24) 



where 



Eq. (|2.15|) is also rewritten 

PuS" = 0. 



(2.25) 



(2.26) 



By this condition, the relation between p^^ and v^^ is de- 
termined as follow : 
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where 



pupa 



^*^pvoifD^ pa- 



(2.27) 



(2.28) 



In this section, we show our numerical results for the 
orbital motion of a test particle, and the wave forms and 
the energy spectra of the emitted gravitational waves. 
To analyze chaotic behavior of a test particle, we use the 
Poincare map and the Lyapunov exponent. To draw the 
Poincare map, we adopt the equatorial plane {9 = 7r/2) 
as a Poncare section and plot the point (r, v^) when 
the particle crosses the Poincare section with < 0. If 
the motion is not chaotic, the plotted points will form a 
closed curve in the two dimensional r-v^ plane, because 
a regular orbit will move on a torus in the phase space 
and the curve is a cross section of the torus. If the or- 
bit is chaotic, some of these tori will be broken and the 
Poincare map does not consist of a set of closed curves 
but the points will be distributed randomly in the allowed 
region. From the distribution of the points, we can judge 
whether or not the motion is chaotic. 



A. Motion of a test particle 



a spinless particle 



As we see in Sec III Bl the equations of the motion of 
a spinless particle are integrable because of the existence 
of the Carter constant. That is chaos does not occur. 
But Johnston found that the orbits of a smnless particle 
can be complicated as if chaos occurred 28] . Integrating 
numerically the equations of motion (|2.5|I - H2.8|I . we show 
a typical complicated orbit in Fig. ^ 

This orbit is called the orbit (a) in this paper. The 
behavior of the orbit (a) looks complicated. However, 
showing the Poincare map in Fig. [21 we find a closed 
curve, which confirms that the orbit (a) is not chaotic. 



2. a spinning particle 

In the previous work [TtL ll^ . the parameter region 
for a spinning test particle in which a particle will show 
chaotic behavior was discussed. We adopt such param- 
eters, e.g. the energy E, the angular momentum J^, 
and the spin parameter 5" of a test particle. Using the 
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FIG. 1: The orbit of a spinless particle with E — 0.968/^, Jz — 
2.0^M,C = 10^2 and a = 1/^2. We choose the ini- 
tial position and velocity as {r,9,(j)) = (lOM, 7r/2,0) and 
(u^tj^u"^) = (0.14,0.03,0.02). This is called the orbit (a) 
in this paper. This orbit looks very complicated as if chaos 
occurred. 




r[M ] 

FIG. 2: The Poincare map of the orbit (a). The plot points 
form a closed curve in the two-dimensional r-v^ plane. We 
also plot the Poincare map of other two orbits with the same 
conserved quantities as those of the orbit (a). This confirms 
that the system is integrable. 




FIG. 3: The orbits of a spinning particle with E = 
0.9328A1, Jz = 2.8/.1M, and S = l.Q^M. We set a = O.SM 
and choose the initial position and velocity as (r, d, (j)) = 
(6.0M,7r/2,0) and u''') = (0.18,0.05,0.07). This is 

called the orbit (b), which is. complicated just as the orbit 
(a). 

0,5 I 1 1 1 1 

(b) 

■ '>!?^. ' ■ "~ ' 

■ :*■-.>:, -i;>s? '7"; 

' " " 

-0.5 I 1 1 1 1 

2 4 6 8 10 12 
r[M] 

FIG. 4: The Poincare maps of the orbit (b) and (c). The plot 
points for the orbit (b) distribute randomly, which confirms 
that it is a chaotic system. For the orbit (c), we will discuss 
it later. 



Bulirsch-Store method ^5Z|, we integrate the equations 
of motion Eqs. (|2.22|) - (|2.24|) and use the constraint equa- 
tions Eqs. lfTTH|l . (|??^ . lfO]|l . (|OT|) to check the accu- 
racy of our numerical integration. We find that the rel- 
ative errors are smaller than 10""'^^ for each constraint. 
We show the orbits of the spinning particle in Fig. O and 
call it the orbit (b). Comparing Fig. ^ with Fig. |31 we 
cannot distinguish two orbits. However, the difference 
between chaotic orbit and nonchaotic one will be appar- 
ent when we draw the Poincare maps. In Fig. ^ we 
show the Poincare map of the orbit (b) . The plot points 
distribute randomly in Fig. 0] We also calculate the Lya- 
punov exponent A to evaluate the strength of chaos. The 
result is depicted in Fig. |31 As shown in Fig. [S] the 
Lyapunov exponent A is positive, which means that the 
orbit (b) is chaotic. The typical time scale for chaos is 
given by A^^ ^ 18M. Since the average of orbital period 
is 74. 3M, the orbit (b) becomes chaotic just after a few 
revolutions around a black hole. 



A[M 
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"T 1 1 r 



4000 8000 12000 16000 20000 
t[M] 

FIG. 5: The maximal Lyapunov exponent A for the orbit (b). 
The typical time scale for chaos is given by A^^ ~ 18M, while 
the average of orbital period is 74. 3M. 
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B. Gravitational wave forms and energy spectra 

Based on the previous calculation of the orbits, we 
show the wave forms and its energy spectra and analyze 
an effect of chaos on the emitted gravitational waves. To 
estimate the gravitational waves, we use the multipole 
expansion of gravitational field |25l l2a | . Here we briefly 
summarize the multipole expansion method. According 
to |2^, the gravitational waves are represented in 
TT gauge as 



1=2 m=~l 



(3.1 



where 



00 

Wi„, Xi^ 

VM^ + I) \ * - sin^ ewi 



Im 







Hm 



^ -Xijn/svae svaOWim | , 

^Anl{l + 1) \ * sin6'Xi„ 



with 



/ 92 d 1 92 



(3.2) 



(3.3) 



(3.4) 



(r, 9) is a position of an observer, dim and fim are an 
odd-parity and even-parity mode, respectively. J'™ is a 
multipole of the mass distribution and S*'™ is that of the 
angular momentum. If a gravitational wave source moves 
slowly, we can write down l'™ and 5'™ as 



rim levr ^ (; + i)(/ + 2) ^^ [ .3 

' =(27TT)!!i 2(^-1); ) y-™^'™^^-' 



(3.5) 



, ^ 32^ / (/ + 2 )(2^ + 1) ^ ^ 
(2Z + 1)!!V2(;- 



{2l + l)\\\2{l-l){l + l), 
X / e,p,a;p(-To,)i;'-^'""V'-id^x, (3.6) 



where Y,' 1^'™* Q^jig^j ^j^e pure orbital harmonic func- 



tion 



I' 1 



I ,lm / 



= E E (i^'™"m'i/m)er'V''"', 



m' — — /' m" — — 1 



(3.7) 
(3.8) 
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(ii) I = 3 
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FIG. 6: The gravitational wave forms of the -I— and the x- 
modes for the orbit (a). The observer is on the equatorial 
plane. Figs, (i) and (ii) show the wave forms of I— 2 and 3, 
respectively. 



Giving the information of a particle orbit, we can eval- 
uate hjj'" for each I. The gravitational waves from a 
chaotic orbit of a spinning particle may contain higher 
multipole moments {I > 2) than quadrupole (I — 2). This 
is the reason why we analyze the gravitational waves not 
only / = 2 but also 1 — 3. Note that because the orbit 
considered here is very rclativistic, the multipole expan- 
sion may not be valid and then it will provide only a 
qualitative feature. 

In Fig. El and Fig. [7| we show the gravitational wave 
forms for the orbits (a) and (b) , respectively. These fig- 
ures reveal two important points. One is that in the wave 
forms, there looks some difference between two orbits (a) 
and (b), but it may be difficult to distinguish which one 
is chaotic. Just as reference, we depict the wave form for 
a circular orbit in Fig. |S1 which is quite regular. The sec- 
ond points is that the amplitudes of the octopole wave 
(I — 3) are smaller than those of the quadrupole wave 
{I ~ 2). Comparing the peak values, we find that the 
ratio of the amplitude of ^ = 3 to that of / = 2 is about 
30 percents. 

In Fig. |51 we show the energy spectrum of the 
quadrupole gravitational wave {I — 2) from the orbit (a). 

In Fig. IIUI we show the energy spectrum of the gravi- 
tational wave(^ = 2, 3) from the spinning particle. From 
Fig. Inland Fig.^| we find that there is a clear difference 
in the energy spectra of the gravitational waves between 
the orbit (a) and (b). 

How about the orbit of a spinning particle which looks 
regular? In order to understand a role of chaotic behav- 
ior, we also analyze the orbit (c), which Poincare map 
is nearly a closed circle (the inner circle map in Fig. 
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(ii) I = 2 
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FIG. 7: The gravitational wave forms of the +- and the x- 
modes for the orbit (b). The observer is on the equatorial 
plane. Figs, (i) and (ii) show the wave forms of I— 2 and 3, 
respectively. 
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FIG. 
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The wave forms for a circular orbit at r = 2.45M. 
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FIG. 10: The energy spectra of the gravitational wave (/ = 
2,3) from the spinning particle. (i)-(ii)The energy spectrum 
of the gravitational wave {I = 2, 3). The values in these case 
are normalized by the peak value of the spectrum / = 2. 
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0.02 0.025 
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FIG. 9: A energy spectrum of the gravitational wave from 
the orbit (a). Each value is normalized by the maximum peak 
value. 



Although this map looks closed in Fig. 01 we find that it 
is not really closed when we enlarge it (see Fig. We 
show the wave forms for the orbit (c) in Fig. E| We also 
present the spectrum in Fig. 1131 

From the figures, we may conclude that the spectrum 
of gravitational waves from nonchaotic orbits contains 
only discrete characteristic frequencies. On the other 
hand, the spectrum for the chaotic orbits contains the 
various frequencies. It seems to be a continuous spec- 
trum with finite widths. The ratio of the octopole wave 
energy for the orbit (b) to that of the quadrupole wave 
is about 12.4 percents. Therefore the higher pole mo- 
ment than the quadrupole moment does not contribute 
so much to the emitted gravitational waves, although the 
orbit is chaotic. 
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FIG. 11: (i) The Poincare map of the orbit (c). We also plot 
the Lyapunov exponent (ii), which value is smaller than that 
of the orbit (b) but is still positive (A ~ 0.02M). 
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FIG. 12: The gravitational wave forms of the +- and the 
X -modes for the orbit (c). The observer is on the equatorial 
plane. Figs, (i) and (ii) show the wave forms of 1= 2 and 3, 
respectively. 
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FIG. 13: A energy spectrum of the gravitational wave from 
the orbit (c). This is similar to the orbit (a). 



IV. SUMMARY AND DISCUSSIONS 

In this paper, we study the gravitational waves emitted 
from a chaotic dynamical system. As a concrete example, 
we analyze the motion of a test particle going around the 
Kerr black hole. We confirm that the orbit of a spinning 
particle can be chaotic. Using the multipole expansion 
method of a gravitational field, we evaluate the gravita- 
tional waves from a chaotic and a non-chaotic orbit in 
order to analyze the effect of chaos on the emitted grav- 
itational waves. As the results, there are not so much 
difference for the gravitational wave forms. However, the 
energy spectra of the gravitational waves show a clear 
difference. For a chaotic orbit, we find a continuous en- 
ergy spectrum with several peaks. While, in the case 
of nonchaotic orbit, the spectrum contains the discrete 
characteristic frequencies. We also find that the higher 
pole moments than quadrupole moment of the system 
do not contribute so much to the emitted gravitational 
waves even for a chaotic orbit. 

When the gravitational waves are detected and the en- 
ergy spectrum is determined by observation, not only the 
astrophysical parameters ,e.g., the mass, the angular mo- 
mentum, and the spin are determined but also some fun- 
damental physics such as relativistic nonlinear dynamics 
could be discussed. 
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As a future work, we have to study the gravitational 
waves in the dynamical system in a full relativistic ap- 
proach. For the first step, using the Newmann-Penrose 
formalism (s^, we extend the present evaluation (the 
multipole expansion) to a full relativistic perturbation 
analysis. 
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